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Quark number susceptibilities from resummed perturbation theory 
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We evaluate the second and fourth order quark number susceptibilities in hot QCD using two 
variations of resummed perturbation theory. On one hand, we carry out a one-loop calculation 
within hard-thermal-loop perturbation theory, and on the other hand perform a resummation of 
the four-loop finite density equation of state derived using a dimensionally reduced effective theory. 
Our results are subsequently compared with recent high precision lattice data, with which they are 
seen to agree down to temperatures very close to the transition region. 
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Introduction. One of the most important challenges in 
the equilibrium thermodynamics of QCD is to develop 
nonperturbative tools to access the region of nonzero 
quark densities, addressing questions such as the exis- 
tence and location of a critical point on the phase dia- 
gram. Barring a solution to the sign problem of lattice 
QCD, the leading method to determine the finite density 
equation of state (EoS), or the behavior of the pressure 
as a function of quark chemical potentials (i, is through 
the evaluation of quark number susceptibilities, 



Xijk...(-L) = 



dfii d(j,j dfik ■ 



(1) 
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where the indices i, j, k, ... refer to different quark flavors. 
These functions carry information about the response of 
the system to nonzero density, yet can be determined 
on the lattice without problems; for examples of recent 
studies, see e.g. [TJ [2] and references therein. The appli- 
cability of these results to the determination of the EoS 
at /i ^ is ultimately restricted only by the convergence 
of the expansion of the pressure in powers of /x/T. 

While a quantitative description of the quark gluon 
plasma near its transition temperature T c clearly neces- 
sitates the use of nonperturbative techniques, it is also 
very interesting to study, to what extent the behavior 
of the quark number susceptibilities can be understood 
using analytic weak coupling methods. First, unlike lat- 
tice calculations, perturbation theory works best at very 
high temperatures and thus offers a way to connect the 
results obtained around T c to arbitrarily high energies. 
More importantly, perturbative calculations are easily 
generalizable to finite density, and are not constrained to 
the region of small n/T. Finally, due to the disappear- 
ance of purely gluonic contributions to the quark number 
susceptibilities, these quantities are expected to display 
improved convergence properties in comparison with the 
pressure itself. Indeed, extensive work on susceptibilities, 
and more generally the chemical potential dependence of 
the pressure, has been carried out within unresummed 
perturbation theory [2111], the hard-thermal-loop approx- 
imation [5H5] , the analytically tractable large- Nf limit of 
QCD [TQim, and even the gauge/gravity duality [12]. 
In addition to this, the applicability of dimensional re- 
duction (DR) to finite densities has been investigated 



in |13j . and the behavior of the susceptibilities deter- 
mined through a nonperturbative DR study in |14] . 

While many of the perturbative calculations listed 
above showed reasonably good agreement with lattice re- 
sults existing at the time of their publication, the numer- 
ically significant corrections present in recent high pre- 
cision lattice data [T] [5] clearly call for a re-examination 
of the issue. On top of this, the past years have wit- 
nessed important progress in the resummation of high- 
temperature perturbation theory on multiple fronts. In 
hard-thermal- loop perturbation theory (HTLpt) [15] . a 
recent evaluation of the partition function of hot QCD 
up to three-loop order has demonstrated dramatically 
improved convergence properties [161 117] , and the agree- 
ment between the HTLpt and lattice results has subse- 
quently been observed to be very good down to 2 — 3 T c 
(for the relevant lattice data, see e.g. [TBH2D])- In addi- 
tion, the same framework has been applied to the case of 
finite density and zero temperature, albeit at lower or- 
ders [21] [22] . At the same time, it was shown in [23] [24] 
that a simple resummation of the soft, three-dimensional 
contributions to the four-loop EoS of hot QCD [55] is 
enough to considerably decrease its renormalization scale 
dependence, resulting in excellent agreement with lattice 
data. It should be very interesting to see, what kind of 
an effect these new techniques have when applied to the 
evaluation of quark number susceptibilities. 

In the present paper, our objective is simple: We want 
to apply state-of-the-art resummation techniques to the 
determination of the second and fourth order quark num- 
ber susceptibilities, and compare the results to the most 
recent lattice data available. To this end, we address two 
separate calculations: First, we employ HTLpt to deter- 
mine the susceptibilities at one-loop order, and then ap- 
ply the resummation scheme of [23J [M] to the four-loop 
|U 7^ EoS derived in [3] to obtain DR results for the 
same quantities. All of our calculations are performed 
within the MS scheme of dimensional regularization, de- 
noting the renormalization scale by A. Furthermore, we 
will in this paper specialize to the phenomenologically 
relevant case of three dynamical quark flavors, which for 
simplicity are all taken to be massless. We have, however, 
explicitly verified that keeping the leading order strange 
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quark mass dependence in the results only affects them 
in any noticeable way at the very lowest temperatures. 

HTL perturbation theory. Hard-thermal-loop pertur- 
bation theory is a reorganization of perturbative expan- 
sions within thermal QCD. The Lagrangian density of 
the theory is written in the form 

^HTLpt = (£qcd + £htl) + A£htl , (2) 

where £qcd is the usual Lagrangian density of the the- 
ory, £htl the HTL improvement term, and 6 a formal ex- 
pansion parameter introduced for bookkeeping purposes. 
The last part of the Lagrangian, A£htl, on the other 
hand contains counterterms, which are necessary to can- 
cel the ultraviolet divergences introduced by the HTLpt 
reorganization. For full QCD with dynamical quarks, the 
(gauge invariant) HTL improvement term reads 



HTL 



S)m D Tr 



F, 



// a 



a n 3 



v y 



-(i 



Nf 

8)iy m 2 



S 



(yD) 
yD 




4>f, (3) 



where D^ L = d^ L — igA^ denotes a covariant derivative (in 
the appropriate representation), y = (l,y) is a light-like 
four vector, () y represents an average over the direction 
of y, and mr> and Tn q t are the Debye mass and fermion 
thermal mass parameters. Note that m q j carries depen- 
dence on the flavor index /, running from 1 to Nf = 3. 

HTLpt is formally defined as an expansion of phys- 
ical quantities in powers of S around <5 = 0, implying 
that already at its leading order 5° one is dealing with 
dressed propagators that incorporate important plasma 
effects, such as Debye screening and Landau damping. 
The starting point of HTLpt is thus an ideal gas of mas- 
sive quasiparticles, which can be identified as one of the 
main reasons for its success. At higher orders, the ex- 
pansion in S generates dressed vertices as well as higher 
order terms that ensure that there is no overcounting of 
Feynman diagrams. 

In practice, physical observables are calculated within 
HTLpt by expanding the corresponding expressions in 
powers of 5, truncating at some specified order, and then 
setting 5 = 1. If it were possible to carry out the expan- 
sion to all orders, the final result would be independent 
of the parameters mu and rn q j. At any finite order in (5, 
some dependence however remains, and a prescription for 
choosing their values is required. Optimally, they should 
be determined via a variational condition for the ther- 
modynamic potential, but at leading order, this would 
lead to mj) = m q.f — 0- I n our calculation, we there- 
fore identify these parameters with their weak coupling 
values, 
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where we have kept the number of colors N c arbitrary. 

After the definitions above, the one-loop HTLpt deter- 
mination of the EoS follows to a large extent the /i = 
work of [15] . In the high temperature limit, we further 
perform an analytic expansion of our result in powers of 
ttid/T and m g /T, resulting in 



d A n 2 T 4 
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where we have denoted d A = N 2 — 1 as well as introduced 
the dimensionless parameters rh = ^ , etc. Results for 
the quark (or baryon) number susceptibilities are then 
obtained by taking derivatives of this expression with re- 
spect to the chemical potentials, and in the end finally 
setting fi = 0. 

Dimensional reduction. To date, the unresummed 
weak coupling expansion of the QCD pressure has been 
worked out up to and partially including its four-loop 
order, both at zero density [2"5Tf2"T] and at /i ^ [3]. A 
very useful tool in these calculations has turned out to be 
the three-dimensional effective theory electrostatic QCD 
(EQCD), the partition function of which contains the 
contributions of the soft and ultrasoft momentum scales 
(gT and g 2 T, respectively) to the corresponding quan- 
tity in the full theory [28, 29 . In practice, one writes the 
pressure of the four-dimensional theory in the form 



PQCD = PHARD + TpEQCD 



(7) 



where phard is defined as the result of a strict loop ex- 
pansion in the full theory, obtained by letting dimen- 
sional regularization regulate both the UV and IR di- 
vergences. At the same time, Peqcd is obtained from 
the partition function of EQCD, which one can evaluate 
using a combination of perturbative [261 !27| and nonper- 
turbative [301 El] tools. 

Formally, EQCD is a three-dimensional Yang-Mills 
theory coupled to an adjoint Higgs field originating 
from the zero Matsubara mode of the four-dimensional 
temporal gauge field. The theory is defined by the La- 
grangian density 

£eqcd = S3 2 j \ Tr[F„f + Tr[(AA,) 2 ] + m| Tr[A%] 

+i( Ti[Al] + A E Tr[^] }+6£ E , (8) 
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where we have assumed iV c = 3 (for larger N c we would 
have two independent quartic terms for the Aq field) , and 
where the last term SCe stands for a series of higher order 
non-renormalizable operators that start to contribute to 
the EoS only beyond 0(g e ). The theory is parametrized 
by four constants: The three-dimensional gauge coupling 
<73, the electric screening mass mj, the cubic coupling 
C ~ E/M/ ( see E2] for details), as well as the quartic 
coupling Ae- All of these parameters have expansions 
in powers of the four-dimensional gauge goupling g, and 
their values have been determined to the accuracy re- 
quired by the four-loop evaluation of the EoS, some even 
beyond this (see e.g. [33]). 

As discussed in [21], the above way of writing the 
full theory pressure suggests a very natural resummation 
scheme: While the unresummed weak coupling expan- 
sion is obtained by expanding the (perturbatively deter- 
mined) EQCD partition function in powers of the four- 
dimensional gauge coupling g 7 one may alternatively sim- 
ply skip this last step and keep peqcd a function of the 
effective theory parameters, writing 



7>eqcd = p M + Pg, 



(9) 



where the functions pm and pc can be read off from 
eqs. (3.9) and (3.12) of 0]. In [23], this procedure was 
observed to lead to a considerable improvement of the 
convergence and renormalization scale dependence of the 
full theory pressure at zero chemical potential. It can, 
however, be applied to the case of the finite density pres- 
sure or the quark number susceptibilities with equal ease, 
which is what we have implemented in our calculations. 
An important step in this in principle straightforward ex- 
ercise is to use the effective theory parameters in a form, 
where they have been analytically expanded in powers of 
fi/T; cf. appendix D of [J and appendix B of [51] . We 
refrain from writing the resulting, very long expressions 
here, but simply display the result of the procedure in 
the plots to follow. 

Choice of parameters. Before proceeding to a quanti- 
tative comparison of our predictions with lattice data, we 
will briefly discuss our choices for the parameters appear- 
ing in the results. These include the values of the renor- 
malization scale A and the QCD scale A^g, in addition 
to which a prescription for determining the form of the 
running gauge coupling must be specified. In both the 
HTLpt and DR calculations, we follow standard choices 
used in the literature, which we summarize below. 

In perturbative calculations of bulk thermodynamic 
observables, the renormalization scale A is typically given 
a value of roughly 2ttT and then varied by a factor of 2 
in order to measure the sensitivity of the result with re- 
spect to this choice. Optimally, the central value should 
result from a presecription such as the Fastest Apparent 
Convergence (FAC) or the Principle of Minimal Sensi- 
tivity (PMS). For the HTLpt result, neither of these is 
however available, and hence the central value is chosen 
as 2irT. In the DR calculation, we on the other hand 
follow a commonly used prescription introduced in [29] 
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FIG. 1. A comparison of our HTLpt (wider, red band) 
and DR (blue band) results for the second order baryon 
number susceptibility xb/T 2 with the lattice results of the 
HotQCD Q] (black dots, extending to T « 250 MeV) and 
Wuppertal-Budapest i2j (WB, green dots) collaborations. 
The bands corresponding to the perturbative results originate 
from varying the values of A and A^jg within the ranges indi- 
cated in the text. Asymptotically, all of the results approach 
the limiting value of 1/3. 



and apply FAC to the three-dimensional gauge coupling 
53, thus obtaining A centra i w 1.445 x 27rT. 

For the dependence of the gauge coupling constant on 
the renormalization scale, we use a one-loop perturba- 
tive expression in the HTLpt result and a two-loop one 
in the DR case. This is in accordance with the usual 
rule that the uncertainties originating from the running 
of the gauge coupling should not exceed those due to the 
perturbative computation itself. Finally, for the choice of 
the QCD scale A^jg- we use a recent lattice determination 
of the value of the strong coupling constant at a refer- 
ence scale of 1.5 GeV [35] • Requiring that our one- and 
two-loop running couplings agree with this, we obtain 
the values of 176 and 283 MeV in these two cases, re- 
spectively. To be conservative, we vary the value of Ajje- 
around these numbers by 30 MeV, which is somewhat 
larger than the reported lattice error bars. 

Results. In Fig. [T] we display our results for the sec- 
ond order baryon number susceptibility \b = d 2 p/dn 2 B , 
which to a very good accuracy satisfies the relation 
Xb = Xmu/3 and for which most of the lattice data has 
been derived. As the widths of the red and blue bands — 
corresponding respectively to the HTLpt and DR results 
- demonstrate, the dependence of our results on the 
renormalization scale and the value of Aj^g is rather mild. 
For instance, a comparison of the DR band with the un- 
resummed four-loop result of [3] shows a reduction of the 
uncertainty by a factor of nearly 10 in this temperature 
range. Our two results are in addition in reasonably good 
agreement with each other, considering that the current 
HTLpt result is only of one-loop order. A comparison 
with the recent continuum extrapolated lattice data of 
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FIG. 2. A comparison of our results for the fourth order 
susceptibility Xmuu with the p4 lattice data of [31] (RBC-B) 
and [37j (Cheng et al.). In the T — > oo limit, this quantity 
approaches the value 6/ty 2 « 0.61. 



both the HotQCD [T] and Wuppertal-Budapest col- 
laborations further shows very good agreement down to 
temperatures of roughly 250 MeV, or less than twice T c . 

Moving on to fourth order susceptibilities, we no longer 
have continuum extrapolated lattice data at our disposal. 
In Fig. [2] we nevertheless compare our results for the light 
quark number susceptibility Xuuuu to the lattice data of 
[5S1 157] . obtained with the p4 action. This time, the 
perturbative results no longer overlap with each other, 
but are nonetheless both in reasonable agreement with 
the lattice. Our hope is that the emergence of continuum 
extrapolated lattice results will eventually help to clarify 
the situation. 

Conclusions. In the present paper, we have applied 



two types of resummed perturbation theory to the deter- 
mination of the second and fourth order quark number 
susceptibilities in thermal QCD. Our main results are 
displayed in Figs. [T] and [2j of which in particular the for- 
mer shows impressive agreement with recent lattice data 
down to temperatures surprisingly close to the transi- 
tion region. Although the state-of-the-art perturbative 
determinations of the zero density EoS compare to the 
lattice results quite favorably, an agreement at this level 
is clearly exceptional (for another very successful case, 
see [35] )■ We interpret this behavior as a reflection of the 
fermionic nature of the susceptibilities; after all, due to 
the chemical potential derivatives appearing in eq. ([I]), 
the most poorly convergent, purely gluonic contributions 
to the EoS altogether drop out from this quantity. 

The most important virtue of weak coupling methods 
is clearly their versatility. Indeed, as soon as the quark 
number susceptibilities for the three-flavor case treated 
here have been computed, it is next to trivial to extend 
the results to other theories of interest, such as two- flavor 
or quenched QCD (or even QCD with a different number 
of colors), as well as to other quantities, such as higher or- 
der susceptibilities or the pressure as a function of chem- 
ical potentials. All of these cases are examples of direc- 
tions we will pursue in a forthcoming publication [39], 
hoping that the results will find phenomenological appli- 
cations in the study of the current and future heavy ion 
data from RHIC, LHC and FAIR. 
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